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The ability to constrain dark energy from the evolution of galaxy cluster counts is limited by the 
imperfect knowledge of cluster redshifts. Ongoing and upcoming surveys will mostly rely on redshifts 
estimated from broad-band photometry (photo-z's). For a Gaussian distribution for the cluster 
photo-z errors and a high cluster yield cosmology defined by the WMAP 1 year results, the photo-z 
bias and scatter needs to be known better than 0.003 and 0.03, respectively, in order not to degrade 
dark energy constrains by more than 10% for a survey with specifications similar to the South 
Pole Telescope. Smaller surveys and cosmologies with lower cluster yields produce weaker photo-z 
requirements, though relative to worse baseline constraints. Comparable photo-z requirements are 
necessary in order to employ self-calibration techniques when solving for dark energy and observable- 
mass parameters simultaneously. On the other hand, self-calibration in combination with external 
mass inferences helps reduce photo-z requirements and provides important consistency checks for 
future cluster surveys. In our fiducial model, training sets with spectroscopic redshifts for ~ 5% — 
15% of the detected clusters are required in order to keep degradations in the dark energy equation 
of state lower than 20%. 



I. INTRODUCTION 

The abundance of clusters of galaxies as a function of 
their mass and redshift is potentially a powerful cosmo- 
logical probe. The sensitivity to the underlying cosmol- 
ogy comes from the dependence of the abundance on the 
comoving volume element and, more importantly, from 
the exponential sensitivity of the cluster mass function 
to the amplitude of linear density perturbations. Both 
of these features depend upon the matter content of the 
universe and the underlying theory of gravity. Possible 
applications include constraints on theories of modified 
gravity [l|, neutrino masses @, the total matter density 
fijvf and the amplitude of linear fluctuations as @, [If. 
Furthermore, cluster counts as a function of redshift offer 
a promising technique to explore and constrain dark en- 
ergy parameters because of the suppression in the growth 
of perturbations during the acceleration epoch. In prac- 
tice these studies are done by comparing observations to 
theoretical predictions from simulations as a function of 
cosmology. 

However, there are many observational challenges to 
the use of clusters to constrain cosmology. First, there 
are different techniques of cluster detection, each one at- 
tempting to obtain cluster samples as complete and clean 
as possible. Typical cluster finding methods explore sig- 
nals such as the Sunyaev-Zel'dovich (SZ) flux decrement, 
X-ray temperature, X-ray surface brightness, overdensi- 
ties in space and color from optical observations and the 
weak lensing shear. When comparing the properties of 
the observed samples to simulation predictions, knowl- 
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edge of the selection function is essential since observa- 
tional effects particular to each cluster finder need to be 
properly accounted for. 

Next, cluster masses must be estimated in ways that 
may or may not be tied to the cluster finder employed. In 
general, the mass is not a direct observable and needs to 
be obtained through the relation between an observable 
proxy and mass, the observable-mass relation. Uncer- 
tainties in mass conversion can lead to degenerate effects 
that destroy most of the information in cluster counts if 
not well calibrated. 

In order to overcome this degeneracy a set of so-called 
self-calibration techniques has been developed recently. 
By requiring consistency between number counts and 
other cluster properties, it is possible to solve simulta- 
neously for cosmology and observable-mass parameters. 
This can be accomplished by follow-up of a small cluster 
sample @, by using information on the clustering prop- 
erties of clusters from their power spectrum [j| or their 
sample covariance from counts in cells 0, 0] , and by using 
information from the shape of the observed mass func- 
tion 0, H|. Another approach uses physical models of 
cluster structure @ which is similar in spirit to imposing 
priors on observable-mass parameters. Combining these 
approaches allows for tests of the assumptions underlying 
the individual approaches. 

Finally, cluster redshifts must be estimated. Whereas 
spectroscopic redshifts are very accurate, when dealing 
with large data sets it becomes impractical to obtain 
spectra for large fractions of objects. Alternatively red- 
shifts can be estimated from broad-band photometry in a 
finite number of filter band-passes (see e.g. and ref- 
erences therein). Photometry can be viewed as a coarse 
spectroscopy that probes the most prominent spectral 
features as they move from their rest positions when the 
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object of interest is redshiftcd. Rcdshifts thus estimated, 
known as photometric redshifts (photo-z's), can be effi- 
ciently calculated for millions of objects and can use not 
only colors but any observable that correlates with red- 
shift [n[. 

Since the interpretation of the counts depends sensi- 
tively on an accurate determination of the redshift dis- 
tribution of the clusters and not on the redshift preci- 
sion for an individual cluster, photo-z's are well suited 
to such studies. Clusters detected by the SZ effect, for 
instance, will typically be followed up optically so that 
photo-z's can be calculated. Optical cluster finding algo- 
rithms can choose to derive photo-z's during the cluster 
finding process itself [H, [l3[ or use externally derived 
photo-z's [HI, EH with different implications for the prop- 
agation of photo-z errors. 

In this paper, we study how the knowledge of the clus- 
ter photo-z error distribution can affect the ability to 
use cluster counts to constrain dark energy. Previous 
works have addressed this question for different cosmo- 
logical probes, including cluster counts [l||, supernova 
[l|| , baryon acoustic oscillations [13, EU and weak lens- 
ing tomography [l9l. [20j . In particular, Huterer et al. [l6| 
studied the effect of systematic shifts in centroids of red- 
shift bins on cosmological constraints, in the context of 
perfect knowledge of cluster masses. Here we generalize 
that analysis by considering the full redshift error distri- 
bution and allowing redshift bias and scatter parameters 
to be arbitrary functions of redshift. We also consider 
the photo-z requirements necessary for self-calibration of 
the observable-mass relation when one is simultaneously 
solving for cluster masses and cosmology. 

We start in |jTT] and EjTTT] describing how redshift errors 
affect cluster number counts and their sample covariancc 
respectively. In ^IVI we describe the fiducial models as- 
sumed and the Fisher matrix formalism, which is em- 
ployed in |JV] to study how redshift errors degrade dark 
energy constraints in various cases of interest. Finally, in 
^ VII we discuss the results and conclude. 



II. NUMBER COUNTS 

For a given cosmology, simulations predict the comov- 
ing number density of dark matter halos as a function of 
mass and redshift. For illustrative purposes, we will iden- 
tify these halos with clusters and employ a fit to simu- 
lations for the halo differential comoving number density 
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where a 2 (M; z) = a 2 R (z), the linear density field variance 
in a region enclosing M = 4irR 3 p m /3 at the mean mat- 
ter density today p m . Because shifts in cluster masses 
can modify cluster counts mimicking a change in cos- 
mology, the exponential sensitivity to a will only be a 
benefit in practice if the observable-mass distribution is 



also well known. Fortunately observations and simula- 
tions suggest observable-mass scaling relations that can 
be parametrized in simple forms and allow for reasonable 
degree of calibration [2^, [H, [24[ ■ Following 0] , we take 
the probability of assigning a mass M ohs to a cluster of 
true mass M to be given by a Gaussian distribution in 
InM 
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For simplicity, we will allow the mass bias In M and 
the variance &\ n M to vary with redshift but not mass. 
The redshift dependent average number density of clus- 
ters within the observable mass range M° bs < M ohs < 
M" 1 ?] is given by 
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where x a = x(M° hs ). Notice that this corresponds to 
the cumulative number density above some sharp mass 
threshold in the limit that of nM — > and — > oo. 
The mean number of clusters is then obtained by inte- 
grating the mean number density in the redshift depen- 
dent comoving volume element d 3 x. Let us use spherical 
coordinates to parametrize the spatial position vector as 
x = (r, 6, <p), where r(z) is the angular diameter distance 
to redshift z and (6, 0) parametrize the solid angle f2 such 
that dil = sin 6d6d<j). Since we will only consider flat cos- 
mologies, r(z) coincides with the comoving distance. The 
volume element is given by 
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where H (z) is the Hubble parameter to redshift z. Red- 
shift uncertainties affect the redshift bin size as well as 
the observed angle, distorting the volume element and 
changing the number counts. We take the probability of 
assigning a photo-z 2:P hot to a cluster of true redshift z 
to be also a Gaussian distribution 
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and the redshift bias z 



bias 



z as (z) and variance a 2 



cr 2 (z) are allowed to be arbitrary functions of redshift. 
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FIG. 1: The redshift distribution of cluster counts per deg 2 above a sharp mass threshold M°£ E = lO 14 - 2 ^" 1 M Q is shown (solid 
line) for the fiducial survey properties defined in £|lVI in the WMAP1 cosmology. The effects of photo-z bias and scatter are shown for 
z bias — az — 0.02 at redshift bins [0.2, 0.3] and [0.4, 0.5] respectively. The dashed lines indicate a perfect redshift selection in the absence 
of photo-z uncertainties. Photo-z errors scatter objects up (light gray) and down (dark gray), changing the observed number counts. In 
the limit that <x 2 — ► 0, a positive redshift bias decreases counts in regions where dm/dz is increasing and vice-versa; the sensitivity of the 
counts is then controlled by —d 2 m/dz 2 at the bin. The effect of scatter in each bin border scales as ~ 0.5d 2 m/dz 2 in this limit; the total 
effect then depends on how d 2 m/dz 2 varies from bin to bin and is proportional to Q.5d s m/dz 3 . 



Given a perfect angular selection characterized by an 
angular top hat window (0), the mean number of 

clusters in a photo-z bin defined by zf hot < z phot < zf^ 
is 



dz phot / d 3 xn a W^(n)p{z phot \z) 



d s xn a Wi(x) 
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where the total window function Wj(x) is given by 

W i (x) = W? h (n)F i (z), (9) 



and 
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The total window function in 



Eq. ([9]) takes into account photo-z error information 
carried by Fi(z). Notice that when z blas = a z = 0, 
the distribution p(z phot \z) = S(z phot - z) and Fi(z) 
becomes a top hat W* h (z), which combined with the 
angular top hat produces the full 3d top hat window 
W* h (x) = Wl h (SY)Wf h (z) of a perfect selection. Using 
Eqs. (0), ©, © and (fTU|) the mean number counts in a 



cell defined by the photo-z bin and a solid angle AO can 
be written as 

m Qil = / dndz^f\n a (z)W^(n)F, t (z) 
J H(z) 

dz ~di~2 t erfc ^^ ~ crfc (^+i)] > ( n ) 

where we used the angular top hat window to perform 
the solid angle integral and defined 



dm a = ^r 2 (z)n a (z) 
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The solid angle AO is related to the angular extension 6 S 
of the cell by 

and the approximation is true in the flat sky regime, valid 
for small angle windows. 

It is instructive to consider the sensitivity of number 
counts to the photo-z parameters in some limits. Let us 
consider the simple case of a sharp mass threshold (i.e. 
drop the a = 1 index) and let us assume that the redshift 
bias and variance are smooth functions of redshift. From 
Eq. pip , the relative sensitivity of the mean counts fhi = 
fna=i,t in a photo-z bin of width Szf = zf,"* — zf ot to 
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whereas the sensitivity to the redshift variance is 
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and the derivatives are evaluated at the photo-z bin cen- 
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zf,° )/2. Note that these derivatives are 
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insensitive to the actual value of a 2 in this limit. Con- 
sequently we will employ the variance a 2 instead of the 
rms a z in the forecasts below. 

In general, the sensitivity of the counts to higher mo- 
ments of the photo-z error distribution will depend on 
higher derivatives of the true redshift distribution and 
will tend to be less important. Fig. Q] shows the redshift 
distribution given by Eq. (p~2|) plotted as a function of 
redshift for the fiducial model defined in § IIVI as well as 
the effects of z bla8 and a z on the observed number counts. 
The distribution is shown only up to z = 0.6 for clarity. 

The result of Eq. (fl"4]) is intuitive from the fact that 
a positive redshift bias decreases the counts in regions 
where the redshift distribution is increasing and vice- 
versa. Consider the effect of the redshift bias on the num- 
ber counts in the redshift bin [z2, 23] = [0.2, 0.3] as shown 
in Fig.[Q We have that dm 2 /dz hms ~ (dm% -Smf )/z bias , 
where Sm% is number of upscattered clusters around zi 
and Smf is the number of downscattered clusters around 
Z3. Approximating these by rectangles, we can combine 

*ft* W M > (16) 



2 dz 

with a similar expression for Smi, and obtain Eq. ((T 



The sensitivity to the redshift variance allows a similar 
interpretation. In Fig. [T] we consider the effect of the 
redshift scatter on the counts in the redshift bin [24, Z5] = 
[0.4, 0.5] . In this case we have dm^j da 2 ~ (Srn^ — Smf + 
(5mg — 8mf)/cj 2 . Approximating these now by triangles, 
we have 
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A similar expression for Smf at Z4 
(Sm^ — Smf) 0.5 dm 



<J Z gives 



2cr z dz 



-0.5 



(17) 



(18) 



(Srrig 



dz 2 

Smf) /a 2 then 



The corresponding result for 
leads to Eq. (IB"]). 

The limits in Eqs. (fl~4"|) and (fT5")) are useful to get qual- 
itative insight into the effect of redshift bias and vari- 
ance when they are small and vary smoothly in redshift. 
Quantitatively, the finite values of these parameters as 
well as more general parametrization of these functions 
may cause the sensitivity of the counts to change from 
these special limits. In £TVl wc explore how the sensitiv- 
ity of the counts to photo-z parameters ultimately affect 
dark energy constraints. 



III. SAMPLE COVARIANCE 

The number counts m Q: ,(x) fluctuate in space tracing 
the linear density fluctuations <5(x) induced by large scale 
structure 



m Q)i (x) = m 0ji [l + b a (z)5(-x)] , 



(19) 



where b a (z) is the average cluster linear bias predicted 
from the distribution in Eq.Q 
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and we take a fit to simulations of [2 
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with a c = 0.75, p c = 0.3, and S c = 1.69. From Eq. (JE 
the counts m a j then possess a sample covariance given 
by M 



{(m a ,i - m at i)(mp,j - mpj)) 
d 3 k 
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W*(k)W 3 (k)P(k) 
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which accounts for the clustering of clusters due to large 
scale structure. Here W^(k) is the Fourier transform of 
the window function and for simplicity we assumed that 
b a (z) does not vary considerably within the photo-z bin 
i and can be approximated by b a (z) w b a (zf) = b a ^. 
In the more general case where b a (z) varies considerably 
within the photo-z bin, Wi (k) above would be the Fourier 
transform of Wi(x)b a (z), and it would be harder to find 
an analytical expression for it. We chose not to consider 
this case since the approximation holds well for the small 
photo-z bin size (6z p = 0.1) considered in our fiducial 
model. 

Let the photo-z bin i be at angular diameter distance 
r.i = r(zf) and have width Srt = T(zf,° t ) — r(zf hot ). 

bias / 



Under the assumption that H(z), z s (z) and a z (z) 
also do not change appreciably inside the bin, so that 
H(z) w H(zf) ee Hi and likewise for z h ' ms (zf) = z bias 
and <r z (zf) ee cr z _i, the window Wj(k) is given by (see 
Appendix fA| 
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Since z^ 13S <C Szf in all practical cases, b a ,i and P(k) 
do not carry any strong dependence on the photo-z bias. 
In this case, Eqs. (|2"2"]) and (|2^|) indicate that the sam- 
ple variance Su does not depend on z° las and is expo- 
nentially sensitive to a 2 i . This will bring interesting ef- 
fects when using sample covariance as a signal for self- 
calibration (see jjVj. From Eq. (|2li|) , we see that in the 
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absence of photo-z errors, the window function would 
suppress modes along the line of sight with wavelengths 
A || < Sr = Sz p /H. The presence of photo-z scatter fur- 
ther suppresses modes Am < ira z /H. For our fiducial 
redshift binning of Sz p = 0.1 we expect significant effects 
to appear then if a z > 5z p /tt ~ 0.03. As we will see in 
ifvT this is roughly the value of scatter uncertainty where 
dark energy degradations start to increase considerably. 

IV. FISHER MATRIX AND 
SELF-CALIBRATION 

Given a parametrized model that predicts the number 
counts and their sample covariance, along with a fidu- 
cial choice for the true values of these parameters, the 
Fisher matrix formalism allows us to study the impact of 
redshift uncertainties on dark energy constraints. 

As discussed below, we will consider counts not only in 
photo-z bins, but also in angular cells and observed mass 
bins. To simplify the notation, from now on i will index 
a generalized pixel of redshift, angle and mass whenever 
these are appropriate and Sij will be the corresponding 
sample covariance. The total covariance is the sample 
covariance plus shot variance 

Gij = S i3 + fhiSij . (24) 

For convenience, we arrange the counts per pixel i into 
a vector m = (mi, . . . , ttln ix ) and correspondingly their 
sample and total covarianccs into matrices S and C. The 
Fisher matrix quantifies the information in the counts on 
a set of parameters p a as 0, H3, [H[ 

F af} = rn^C^m p + ^[C^S.aCr^] . (25) 

The first term represents the information from the mean 
counts and the second term carries the information from 
the sample covariance of the counts. When using only 
counts information, we will be using only the first term 
in the Fisher matrix definition, while when using the 
sample covariance to employ self-calibration we will use 
both terms. The Fisher matrix approximates the co- 
variance matrix of the parameters C a p ss [F^ 1 ]^ such 
that the marginalized error on a single parameter is 
a{p a ) = [F _1 ]ia . When considering prior information 
on parameters of a given a(p a ) we add to the Fisher ma- 
trix a contribution of o~ 2 (p a )$af3 before inversion. 

Next we define our fiducial choices for survey prop- 
erties and values for cosmological parameters as well as 
nuisance parameters describing the observable-mass re- 
lation and photo-z's. 

A. Fiducial Model 

We will take a fiducial cluster survey with specifica- 
tions similar to the South Pole Telescope (SPT) Survey 



[29l | : an area of 4000 deg 2 and a sensitivity corresponding 
to a constant mass threshold M°^ s = lO 14 - 2 /i _1 M . We 
divide the number counts into photo-z bins of Szf = 0.1 
out to z%~£ = 2. When using self-calibration from sam- 
ple variance information, we further divide these counts 
into 400 angular cells of 10 deg 2 . Finally, when cm- 
ploying self-calibration from shape of the observed mass 
function, we additionally divide the counts in 5 bins of 
Alog 10 Af obs = 0.2. In SjVD]wc will consider some vari- 
ations on these fiducial choices. 

The obscrvable-mass relation will be parametrized by 
the mass bias M blas and variance of nM as a function 
of redshift. Because the evolution in cluster parameters 
is expected to be smooth in redshift, we will assume in 
our fiducial model that the mass bias has a power law 
evolution given by 

In M h ' ms (z) = A b + n b ln(l + z) (26) 

with fiducial values A b = rif, = 0. For the mass variance 
of M wc assume a Taylor expansion around z = 

a=0 

with fiducial values of nM |fid = (0.25) 2 and B a = 0. 
We have checked that using a different fiducial value of 
af n M |fid = (0.05) 2 does not change the results. For the 
fiducial model we will take a cubic (N a = 4) evolution 
of the mass variance, so our fiducial model will be com- 
prised of 6 nuisance parameters for the observable-mass 
distribution. In WDI we also consider a more general 
parametrization with one value of Af blas and cr 2 nAf in 
each photo-z bin. 

For the photo-z model wc will take the values of z blas (z) 
and <J z (z) in N z different rcdshifts as our nuisance photo- 
z parameters. Values of photo-z parameters at arbitrary 
values of z are then obtained by cubic spline interpo- 
lation. Current photo-z methods can reach accuracies 
of a z ~ 0.1 for field galaxies up to moderate redshifts 
(z ~ 1.4) when using optical filter band-passes [l(| (see 
also [ll[ for a case with lower redshifts). If complemented 
by near-infrared filters, the accuracy can be improved 
and well controlled up to high redshifts (z ~ 2.0). Photo- 
z's of red early-type galaxies are typically even better 
because of their very prominent 4000A break. Finally, 
cluster photo-z's are expected to be further improved 
by averaging the photo-z's of individual cluster galax- 
ies. However, this will only be the case if the presence 
of interlopers is small and well understood. Wc expect 
therefore not only the photo-z errors to be smaller than 
the corresponding errors for field galaxies, but the error 
distribution to be more well behaved. For instance the 
Gaussian assumption of Eq. © should be a good approx- 
imation for clusters with well understood selections. For 
our fiducial photo-z model, we will take N z = 20, i.e. 
one value of z blas and a 2 for each bin Sz = 0.1 of true 
redshift. We will assume fiducial values z blas (;z)|fid = 
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and 



<r g (z)\ ai = 0.03(1 + z) 



(28) 



to account for the fact that the photo-z scatter is ex- 
pected to increase at high redshifts. In WDI . we consider 
a different photo-z model with a constant <7 2 (z)|fid = 
0.02, and show that the results are very similar to the 
ones in the fiducial model. That the results are insensi- 
tive to the fiducial values assumed in a given parametriza- 
tion is a check of the validity of our Fisher matrix ap- 
proach. Obviously, changes in the parametrization and 
in the number of nuisance parameters affect the results. 
The choice of N z = 20, corresponding to one value of 
z blas and a 2 per photo-z bin used to count clusters (recall 
Szf = 0.1) is, although natural, arbitrary. In principle, 
we could choose N z = 10 or 40 and keep the same photo- 
z binning or alternatively use coarser/finer photo-z bins. 
The photo-z bin size is ultimately dictated by the photo-z 
precision that can be achieved, e.g. the photo-z bin size 
should be larger than z blas or a z . The value of Szf = 0.1 
chosen reflects conservative expectations for the photo-z 
precision of future cluster surveys. In WDI wc consider 
alternative models with N z — 10 and Szf = 0.2. 

For the fiducial cosmological parameters, we will as- 
sume a flat universe and take parameter values based 
on the results of the Wilkinson Microwave Anisotropy 
Probe first year data release (WMAP1). In flVDl we will 
also consider a fiducial cosmology based on the third year 
data release (WMAP3). We chose the WMAP1 as the 
fiducial case since this cosmology will place the stronger 
requirements on the photo-z errors. For the WMAP1 
case, the cosmological parameters are the normalization 
of the initial curvature spectrum S^(= 5.07 x 10~ 5 ) at 
k = 0.05 Mpc -1 (corresponding to <j s = 0.91, see j30l|). 
its tilt n{~ 1), the baryon density Qbh 2 (= 0.024), the 
dark matter density fl m h 2 (= 0.14), and the two dark 
energy parameters of interest: its density fiDE(= 0.73) 
relative to critical and equation of state w(= —1) which 
we assume to be constant. The fiducial values are given 
in parentheses. For the WMAP3 cosmology we take 
5{(= 4.53 x 10~ 5 ) at k = 0.05 Mpc -1 (corresponding to 
cr 8 = 0.76), n(= 0.958), n b h 2 (= 0.0223), n m h 2 (= 0.128), 
^de(= 0.76) and w(= —1). We will assume 1% priors 
on all cosmological parameters except for the dark energy 
parameters, which will vary freely. 



V. RESULTS 

Now we apply the formalism developed in i flVl to study 
how dark energy constraints are affected by redshift un- 
certainties. We first consider the case where masses 
are perfectly known and only information from counts 
is used to constrain dark energy. We then progress to 
a more pessimistic case where masses are unknown and 
self-calibration is employed in a joint fit of dark energy 
and obscrvablc-mass parameters. We then consider the 
case where masses arc known, and self-calibration is used 



to reduce the redshift requirements. Finally, we study 
how stable our results are to changes in our fiducial 
parametrization and estimate the requirements on the 
size of calibration training sets necessary in order not to 
degrade dark energy constraints appreciably. 

Given the values of the prior uncertainties in the red- 
shift bias cr(z blas ) and variance <j(<t 2 ), which produce 
the constraint a(w) = a [w\a(z hlas ),a(a 2 )~\ on the dark 
energy equation of state u>, we define the w degrada- 
tion d w relative to a reference case with w constraint 
c(w)| re f < a(w) as 

d w {a{z ),a(a z )) = — 1 (29) 

O"(w)|rof 

and likewise for the Ode degradation df2 DE • Unless oth- 
erwise specified, we will take the reference case to be the 
baseline case of perfect redshifts cr(w)| rc f = a [w\0, 0]. In 
some cases though we will also use this degradation def- 
inition to generally compare two different scenarios such 
as the WMAP1 versus WMAP3 cosmologies and imple- 
mentations assuming masses to be perfectly known versus 
self-calibrated. 



A. Perfect Mass Calibration 

In the case of perfect knowledge of the observable- 
mass relation, the dark energy constraints from clus- 
ter counts alone are quite strong, as long as we have 
a reasonable knowledge of redshifts parameters. In the 
case of perfect redshifts, the baseline constraints are 
<j(w,tt DE ) = (0.033,0.0081) in the fiducial model. 

These constraints are degraded as the uncertainties in 
the redshift bias and variance increase. With counts in- 
formation only, interesting dark energy constraints can 
only be extracted if we have some prior knowledge of both 
redshift bias and scatter uncertainties, i.e. no constraints 
can be extracted if either a(z bms ) — > oo or <t(<j z ) —> oo. 

The top panels of Fig. [2] show the degradations in 
dark energy for finite values of redshift uncertainties. 
The top left panel shows d w and df2 DE as a function of 
a(z hlas ) (with a(a 2 ) = 0) and as a function of a(a z ) (with 
cr(2: blas ) = 0). The top right panel shows contours of 
fixed d w values in the plane defined by values of <r(,z blas ) 
and ^<j(a1). Table U shows constraints and degradations 
in the fiducial model for the case of perfect masses and 
for the self-calibration cases of £|VBI and £|V CI As ex- 
pected, the dark energy constraints are more sensitive 
to the photo-z bias compared to the photo-z scatter by 
about an order of magnitude. 



B. Mass Self-Calibration 

Now we consider the case where we do not know the 
observable-mass parameters and solve for them along 
with dark energy parameters simultaneously by means 
of self-calibration with clustering and shape information. 
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FIG. 2: Dark energy sensitivity to prior knowledge of redshift parameters in the WMAP1 cosmology and fiducial model described in 
ijlVI The top row displays results for the case of perfect masses employing only counts information. The bottom row is for the case 
of marginalized obscrvablc-mass parameters and self-calibration with clustering and shape information. The left panels show percent 
degradation in dark energy parameter constraints with respect to case of perfect redshift knowledge as a function of prior uncertainty in 
z s for fixed o^ and as a function of prior uncertainty in o\ for fixed z blas . The right panels display the respective contours of fixed 
degradation d w in the <r(z blas ) versus \J a(a'^ ) plane. Compared to perfect masses, in the case of self-calibration dark energy degradations 
are weaker due to (i) worse baseline constraints and (ii) decreased degeneracies between dark energy and photo-z parameters, partially 
broken by self-calibration. 



Let us consider first the case of perfect redshifts. Self- 
calibration results are shown in Fig. [3] for our fiducial 
model. For reference, the case of perfect masses is dis- 
played in the black solid ellipsis corresponding to the 
tight baseline constraints of §V Al If no self-calibration 
of the observable-mass relation is used in a joint fit of 



mass and dark energy parameters, no interesting con- 
straints can be extracted. Keeping redshifts fixed, assum- 
ing the functional forms of Eqs. (j2l))) . ((77)) for the mass 
bias and variance and employing self-calibration from 
clustering information in the sample covariance of counts 
produces constraints a(w, fi DE ) = (0.23,0.22). Further 
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Observable-Mass parameters 


Information 


Constraints 


Degradations 


lnM bias 






a(w) 


o-(Qde) 


d w 


do DE 


Known 


Known 


Counts 


0.033 


0.0081 


21% 


17% 


Known 


Marginalized 


Counts 


0.19 


0.21 


90% 


108% 


Marginalized 


Marginalized 


Counts 










Known 


Known 


Counts + Self-Calibration 


0.031 


0.0078 


16% 


15% 


Known 


Marginalized 


Counts + Self-Calibration 


0.072 


0.019 


5% 


3% 


Marginalized 


Marginalized 


Counts + Self-Calibration 


0.12 


0.029 


11% 


8.3% 



TABLE I: Dark energy baseline constraints and degradations for the fiducial model with observable-mass parameters known/marginalized 
and with/without self-calibration. Baseline constraints cr(w) and <t(Qde) are for (x(z blas ) = o{a z z ) = and degradations d w and cin DE 
relative to the baseline are for (r(z blas ) = 0.003 and \J o"(<r|) = 0.03. 



-0.7 




0.6 0.7 0.8 0.9 



FIG. 3: Dark energy constraints for the fiducial model in the 
WMAP1 cosmology. From inner to outer ellipses, the 68% CL re- 
gions arc shown for the case where rcdshift and mass parameters are 
perfectly known (black); the case with self-calibration from sam- 
ple covariance and shape information but fixed photo-z parameters 
(light gray); the case with self-calibration and prior uncertainties 
on photo-z parameters of <r(z blas ) = 0.003 and \/o-(cr'j) = 0.03 
(dark gray). The latter case corresponds to a degradation of 
d(w, Ode) = (11%, 8.3%) with respect to the case of fixed red- 
shifts. 



adding shape information from multiple mass bins al- 
lows a good degree of self-calibration as shown in the 
light gray ellipsis of Fig. [31 the constraints in this case 
are <j(w, Ode) = (0.12,0.029), i.e. there is a degradation 
of d(w, SIde) = (280%, 260%) with respect to the case of 
perfect mass in j]V Al 

As the uncertainty in rcdshift parameters increases, the 
ability to self-calibrate masses is reduced. The dark gray 
ellipsis in Fig. [3] shows how the self-calibrated case above 



is degraded if we take the photo-z parameters to have 
uncertainties of cr(z bias ) = 0.003 and ^a{a 2 z ) = 0.03 si- 
multaneously as opposed to being fixed; in this case the 
constraints are a(w, ^de) = (0.14,0.032), i.e. a degra- 
dation of d(w,il£,E) = (11%, 8.3%). The dark energy 
degradation in the case of self-calibration is shown in the 
bottom panels of Fig. [5] and in the bottom rows of Ta- 
ble [T] Notice that because the baseline constraints (per- 
fect rcdshifts) arc worse with self-calibration compared to 
fixed masses, the degradation relative to the baseline is 
correspondingly weaker. The weaker degradation is par- 
tially because the baseline constraints are worse, but also 
because self-calibration helps break not only the degen- 
eracy between dark energy and mass, but also redshifts. 
In CI we investigate that further. 



C. Redshift Self-Calibration 

If external mass calibrations can partially or totally 
constrain the mass bias or scatter or both, self-calibration 
is still useful to monitor rcdshifts, therefore decreas- 
ing the rcdshift requirements. More generally, self- 
calibration provides important consistency checks that 
theoretical assumptions and observations agree, even 
more so when we have stringent external mass priors. 
External mass calibrations can come for instance from 
follow-up of a small sample of clusters with very well 
measured masses, weak lensing calibrations, simulations, 
etc. 

Let us consider the limiting case where external cali- 
brations put strong priors on the mass bias, but no prior 
on the mass scatter, i.e. cr(lnM blas ) = and <j(crf nM ) = 
oo. In this case, without self-calibration, the constraints 
for perfect redshifts arc <j{w,Qve) = (0.19,0.21) and 
they degrade to a(w,^r>E) = (0.35,0.44), i.e. a degra- 
dation of d(w, JIde) = (90%, 108%) if we impose rcdshift 
priors of cr(z bias ) = 0.003 and y/a{a 2 z ) = 0.03. Employ- 
ing self-calibration the constraints for perfect redshifts 
a(w,n DE ) = (0.072,0.019) degrade to a(w,il DE ) = 
(0.076,0.019) with the same redshift priors, i.e. a degra- 
dation of only d(w, ^de) = (5%, 3%). 



9 



Consider now the limiting case where external priors 
can further constrain the mass scatter to a(af nM ) = 
in addition to the mass bias. Without self- 
calibration, the constraints for perfect redshifts are 
a(w,fl DE ) = (0.033,0.0081) and they degrade to 
(t(w,Qf>e) = (0.039,0.095) if we impose redshift pri- 
ors of cr(z bias ) = 0.003 and ^/ajafj = 0.03. , i.e. 
a degradation of d(w,QuE) = (21%, 17%). Employ- 
ing self-calibration the constraints for perfect redshifts 
<t(w,S] D e) = (0.031,0.0078) degrade to a(w,fl DE ) = 
(0.036, 0.0090) with the same redshift priors, i.e. a degra- 
dation of d(w, ^de) = (16%, 15%). These results are also 
displayed in Table [U Interestingly, in the limit where 
we have no redshift knowledge (er(z blas ) = <j(cr 2 z ) = oo), 
without self-calibration no interesting dark energy con- 
straints are possible, but employing full self-calibration 
we can still obtain a(w, S^de) = (0.11,0.11). 

These results show that self-calibration is very effec- 
tive not only in allowing a joint fit of mass and cosmo- 
logical parameters, but also in reducing the redshift re- 
quirements in those cases. As external priors put con- 
straints on some of the observable-mass relation param- 
eters, but not all of them, self-calibration becomes even 
more important. In the limit where external priors are 
even stronger and masses are perfectly determined, self- 
calibration becomes less important, but it still helps mon- 
itor redshifts. However, even in this case prior redshift 
knowledge is still important in order to extract the survey 
full constraining power. Since any given survey is likely 
to be in between these limit cases, both self-calibration 
techniques and good knowledge of redshifts will be com- 
parably important. 



D. Fiducial Model Dependence 

In this section we consider deviations from our fiducial 
model and their effect on the dark energy degradation 
results. Even though changes in the fiducial values of 
the nuisance parameters describing the observable-mass 
relation and photo-z's have a very small effect on the con- 
straints, changes in their parametrization will obviously 
affect the results. For instance, increasing the number of 
parameters in the mass-observable relation to allow for 
more general models produces worse baseline constraints 
and correspondingly weaker degradations due to photo-z 
uncertainties. The constraints obviously also change if 
the fiducial cosmology and survey specifications deviate 
considerably from the ones assumed in the fiducial model. 
We will first consider a fiducial cosmology based on cos- 
mological parameters from WMAP3 and also changes in 
our fiducial survey properties. Finally we will consider 
deviations in the photo-z model assumed. 



1 1 1 1 1 1 1 1 1 1 
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( S ii /m 1 )«ap3 


- 


( S ii /m i )w«api _ 


i i i i 


i i i i 1 i i i i 1 i i i 



0.5 1 1.5 2 

z 

FIG. 4: The lower value of erg in the WMAP3 cosmology decreases 
the number counts with respect to WMAP1, specially at high red- 
shift. Similarly, the relative importance of sample variance errors 
with respect to the increased shot noise decreases, lowering the 
power of using sample variance as a signal for self-calibration. 



1. Fiducial Cosmology and Survey Properties 

Because of the exponential sensitivity of the cluster 
mass function to the growth of structure, the decrease 
in the value of as for WMAP3 causes the overall num- 
ber counts to be reduced by a factor of ~ 4 with respect 
to the WMAP1 case. Likewise, the reduction in counts 
decreases the relative importance of sample variance er- 
rors with respect to Poisson shot noise. In Fig. 2] we 
show the ratio of the average number counts rrii{z) in 
the WMAP3 and WMAP1 cosmologies, as well as the 
ratio of the corresponding sample variance Su relative to 
Poisson variance rrii. 

Let us first consider the case of perfect knowledge 
of the observable-mass parameters. In the case of per- 
fect redshifts, the baseline constraints for the WMAP3 
cosmology with count information become cr(w, Ode) = 
(0.047,0.010). Even though the WMAP3 total num- 
ber counts are reduced by a factor of ~ 4 compared 
to WMAP1, the baseline constraints degrade only by 
d(w, Ode) = (44%, 25%). If all clusters carried the same 
cosmological information, we would naively expect a fac- 
tor of ~ 2 degradation from the increased Poisson errors. 
However, the rare high mass/redshift clusters are more 
sensitive to cosmology. The fact that the relative reduc- 
tion on the number of these rare clusters is more pro- 
nounced provides extra cosmology information that par- 
tially balances the degradation from the overall counts 
reduction. 
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Case 


Number Counts 


Observable-Mass 


Information 


Constraints 


Degradations 




per deg 2 


parameters 




\ i 


o"(Ode) 




do n[ , 


Fiducial model 


8.0 


Known 


Counts 


0.033 


0.0081 


21% 


17% 


WMAP3 


2.1 


Known 


Counts 


0.047 


0.010 


5.7% 


6.1% 


A/ t ° h bs = W 14 - A3B h- 1 M e 


2.1 


Known 


Counts 


0.050 


0.012 


5.6% 


5.8% 


Area = 1000 deg 2 


2.0 


Known 


Counts 


0.054 


0.014 


4.9% 


3.5% 


phot _ q r 


2.3 


Known 


Counts 


0.063 


0.015 


410% 


780% 


Fiducial model 


8.0 


Marginalized 


Counts + Self-Calibration 


0.12 


0.029 


11% 


8.3% 


WMAP3 


2.1 


Marginalized 


Counts + Self-Calibration 


0.24 


0.048 


7.0% 


6.2% 


M t ° h bs = W 14A3B h- 1 M e 


2.1 


Marginalized 


Counts + Self-Calibration 


0.24 


0.052 


7.3% 


5.4% 


Area = 1000 deg 2 


2.0 


Marginalized 


Counts + Self-Calibration 


0.24 


0.052 


4.0% 


2.9% 


phot _ n c 


2.3 


Marginalized 


Counts + Self-Calibration 


0.28 


0.034 


35% 


42% 



TABLE II: Dark energy baseline constraints and degradations for the fiducial model (WMAP1, A/^J 38 = lO 14 ' 2 /^ 1 Af , Area = 4000 

deg 2 , Zmax =2.0) and various deviations from the fiducial survey properties and cosmology. In all cases the total number of clusters is 
reduced by a factor of ~ 4 with respect to the fiducial case. Baseline constraints are for <r(z blas ) = cr((X 2 ) = and relative degradations 
are for cr(,2 bias ) = 0.003 and y V(cr 2 ) = 0.03. 



Let us now consider the dark energy degradations due 
to photo-z uncertainties. The top panel of Fig. [5] shows 
the degradation of dark energy parameters as a function 
of photo-z uncertainties for the WMAP3 cosmology. No- 
tice that dark energy degradations relative to the baseline 
cases are slightly weaker for the WMAP3 case, since its 
baseline is worse, confirming the trend that surveys with 
smaller yields require less knowledge of redshift param- 
eters for fixed baseline degradations. In other words, in 
the limit of perfect masses for the WMAP3 case, dark en- 
ergy constraints are less sensitive to photo-z uncertainties 
because the best achievable results are worse. 

Now we consider the case where we marginalize 
over the observable-mass parameters and employ self- 
calibration. In this case, with perfect redshifts the con- 
straints for WMAP3 are a(w,fl DE ) = (0.24,0.048), i.e. 
a degradation of cI(w,Qde) = (410%, 370%) with re- 
spect to the case of perfect masses and d(w,Qr>E) = 
(96%, 64%) with respect to the corresponding self- 
calibrated case in the WMAP1 cosmology. Self- 
calibration is based on the ability to divide the cluster 
sample in space and mass. The decrease in the number 
counts for WMAP3 reduces the relative importance of 
sample variance over the increased shot noise, and re- 
duces the power of self-calibration by clustering. In ad- 
dition, the ability to split the sample in mass bins is de- 
creased because we run out of clusters at high masses, and 
self-calibration from shape information is also reduced. 

The bottom panel of Fig. \5\ shows the dark en- 
ergy degradation as a function of photo-z uncertainties 
for the case of self-calibration in the WMAP3 cosmol- 
ogy. Relative to the case of perfect masses and to the 
self-calibrated case of WMAP1, most degradations for 
WMAP3 are weaker, as expected from the worse base- 
line constraints. However the ct(Ode) degradation as a 
function */ ct(ct|) is stronger relative to both cases. That 
happens because self-calibration uses sample covariancc 



information, which is exponentially sensitive to . With 
the reduced relative importance of sample variance in 
WMAP3, changes in are now relatively more degen- 
erate with cosmology and observable-mass parameters. 

To make the point clearer, let us consider the 
case where we employ self-calibration only from shape 
information, but not clustering. In this case the 
self-calibration constraints for perfect redshifts are 
<t(u>,Ode) = (0.26,0.11), i.e. the w constraint is nearly 
identical to the full self-calibration case, but the Ode 
constraint is still ~ 2 times higher. Keeping the photo- 
z bias fixed cr(z blas ) = and applying a prior on the 
photo-z variance of \J c(er^) = 0.2, the constraints de- 
grade to cr(w, Ode) = (0.87,0.17), i.e. a degradation 
of d(w, Ode) = (240%, 54%). With full self-calibration 
the corresponding degradation is (30%, 110%). Therefore 
whereas clustering self-calibration from sample variance 
does not help the w constraint to further decrease con- 
siderably from its value with shape self-calibration only, 
it does make this constraint less sensitive to photo-z un- 
certainties. On the other hand, clustering self-calibration 
improves the Ode constraint, but leaves it more unsta- 
ble to photo-z scatter uncertainties. These features tell 
us that, when employing self-calibration, knowledge of 
the photo-z scatter is essential, specially in a situation of 
fewer clusters such as in the WMAP3 cosmology. 

For completeness, we consider other survey modifica- 
tions that lead to roughly the same reduction in the to- 
tal number counts as the WMAP3 cosmology relative to 
WMAP1. Table ILT1 shows results for the fiducial model 
along with the WMAP3 case and changes in the mass 
threshold M°^ s , total survey area, and maximum survey 
redshift z^°^. All cases lead to a reduction in the total 
counts by a factor of ~ 4 relative to the fiducial model. 
In the case where we change the maximum redshift to 
z max = 0-5, we also reduce the number of photo-z pa- 
rameters to N z = 5 so as to keep one value of z hlas and 
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FIG. 5: Degradation in dark energy parameters similar to Fig. [2] 
but for the WMAP3 cosmology. 



a\ per photo-z bin of Szf =0.1. Except for the change in 
z maxj which makes dark energy constraints and degrada- 
tions much worse, all other cases produce results roughly 
similar. Obviously although these changes produce sim- 
ilar counts reduction, they all remove different clusters. 
For instance, the WMAP3 reduction is more pronounced 
at high rcdshifts, whereas the change in M°^ s removes all 
low mass clusters and the change in survey area reduces 
the numbers of all clusters equally. 

In the case of perfect masses and counts only, except 



for the z ma x reduction case, the cases with worse baseline 
constraints have weaker degradations; the more drastic 
change of reducing zf~j£ has a much higher impact be- 
cause with the total removal of high redshift clusters, the 
leverage to probe the evolution of the growth of structure 
is severely diminished. 

In the case of self-calibration, many factors affect the 
degradations. For the fiducial model and the cases of 
decrease in area and zg^f, the degradations are smaller 
compared to the case of only counts, because their base- 
line constraints are much worse and self-calibration is 
fairly efficient. Even though the degradation in the case 
of 2:^°* reduction is still the highest, it is many times 
smaller than in the case of counts only, because with this 
small redshift range, self-calibration is able to break most 
of the high degeneracy between mass, redshift and cos- 
mology. For the cases of WMAP3 and reduction in M t ° bs , 
the degradations are larger compared to counts only and 
perfect masses. In the WMAP3 case self-calibration is 
harder because of the reduction in erg, which lowers the 
value of Su. In the case of reduction in M°^ s , the removal 
of the low mass clusters makes self-calibration from shape 
information harder. In any case, we confirm the trend 
that more ambitious surveys will be able to extract more 
information from self-calibration techniques but will also 
require more redshift knowledge for a full realization of 
their constraining power. 



2. Fiducial Observable-Mass Parameters 

Now wc return to the WMAP1 cosmology case, but 
we change our fiducial paramctrization of the observable- 
mass relation. We allow the mass bias and variance to 
be general functions of redshift, i.e. instead of assum- 
ing the functional forms in Eqs. (|26|) and (|27|) . we now 
have one value of A/ bias and 

"lnM m 6ach photo-z bin 
of Szf = 0.1. For the case of fixed mass parameters, 
the dark energy constraints obviously do not change, but 
they become a(w,Q^,E) — (0.22,0.030), i.e. a degrada- 
tion of d(w,f2DE) = (80%, 4%) with respect to the fidu- 
cial model with self-calibration. 

In Fig. [6] we show the corresponding dark energy degra- 
dation from the baseline constraints as a function of 
photo-z uncertainties in this case. Now a{w) has a worse 
baseline, and its degradation is less sensitive to er(z blas ) 
but more sensitive to a(a^) compared to the fiducial 
model. On the other hand, since the baseline ct^de) ba- 
sically did not change, meaning that it is making strong 
use of self-calibration, its degradation as a function of 
cr(z blas ) is nearly identical, but is much stronger as a func- 
tion of crial). In this conservative case of 40 parameters 
for the observable-mass relation, where self-calibration is 
being used in its full power, if we are interested in mod- 
erate dark energy degradations (~ 10% — 20%), both 
cr(z blas ) and <j(ct^) are important. As we transit to the 
case with fewer parameters, cr(z blas ) becomes relatively 
more important. 
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FIG. 6: Degradation in dark energy parameters for the case where 
M blas and af n M are free functions of redshift as opposed to hav- 
ing the functional forms of Eqs. (1261) . l|27|l assumed in our fiducial 
model. 



3. Fiducial Photo-z Parameters 

Finally, we consider the case where Eq. (|28|) in our 
fiducial model is replaced by a fiducial photo-z scatter 
of cr z |fid = 0.02 constant in redshift. Let us consider 
the case of perfect masses first. In the case of per- 
fect rcdshifts the baseline constraints in this case are 
cr(w, D e) = (0.032,0.080), i.e. they remain nearly the 
same as in our fiducial case (~ 1% improvement). In the 
top panel of Fig. [7] we show dark energy degradations as 
we allow for uncertainties in photo-z parameters. Since 
the z blas model is still the same, degradations as a func- 
tion of a(z hlas ) are identical to the fiducial case. As a 
function of ^<j(<t z ) they are only slightly weaker. 

In the case of self-calibration, the baseline constraints 
are <r(w, ^de) = (0.11,0.026), i.e. they improved 
by (6%, 10%) with respect to the corresponding self- 
calibrated case in the fiducial model. Now the improve- 
ment is more noticeable than in the case of perfect masses 
since we use sample variance to self-calibrate mass pa- 
rameters. In the bottom panel of Fig. [7] we show dark 
energy degradations in this case. Again, degradations as 
functions of cr(z blas ) are identical to the fiducial model. 
As a function of \fo\p^) degradations become stronger 
at lower values of <j(u z ), since we are now in a better 
baseline. At higher values of <j(ct z ), degradations become 
weaker, since self-calibration is now more efficient. 

Finally, let us consider the effect of changing the num- 
ber of photo-z parameters N z and/or the photo-z bin 
size Sz p . In Table IIIII we show the fiducial model and 
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FIG. 7: Dark energy parameters degradation for a fiducial photo- 
z scatter of o- z |fj d = 0.02 constant in redshift as opposed to the 
fiducial model of Eq. J28}. 



variations in the photo-z model that include changes in 
the fiducial values of the photo-z scatter, N z and 5z p . 
The degradations arc shown for priors of <r(z blas ) = 
0.003^/^/20 and o(pl) = (0.03) V^/ 20 - The scal- 
ing of the priors with N z comes from Eqs. (f3"U)) . (f3~Tj) and 
the fact that the number N spcc of spectroscopic calibra- 
tors in a given redshift region is inversely proportional to 
the number of redshift bins iV spec oc 1/N Z . These scaled 
priors then reflect constraints per fixed bin Sz = 0.1 of 
true redshift independently of the value of N z used. No- 
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Case 




Observable-Mass 


Information 


Constraints 


Degradations 








parameters 




a(w) 


ct^de) 




flf) n[ , 
1J DE 




Fiducial model 




Known 


Counts 


0.033 


0.0081 


21% 


17% 




cr z | fid = 0.02 




Known 


Counts 


0.032 


0.0080 


12% 


13% 


N z 


= 10 and 8z p = 


0.1 


Known 


Counts 


0.033 


0.0081 


15% 


15% 


N z 


= 20 and 8z p = 


0.2 


Known 


Counts 


0.033 


0.0083 


51% 


52% 


N z 


= 10 and Sz p = 


0.2 


Known 


Counts 


0.033 


0.0083 


14% 


14% 




Fiducial model 




Marginalized 


Counts + Self-Calibration 


0.12 


0.029 


11% 


8.3% 




a z | fld =0.02 




Marginalized 


Counts + Self-Calibration 


0.11 


0.026 


15% 


30% 


N z 


= 10 and Sz p = 


0.1 


Marginalized 


Counts + Self-Calibration 


0.12 


0.029 


13% 


14% 


N z 


= 20 and Sz p = 


0.2 


Marginalized 


Counts + Self-Calibration 


0.13 


0.031 


22% 


6.1% 


N z 


= 10 and 5z p = 


0.2 


Marginalized 


Counts + Self-Calibration 


0.13 


0.031 


11% 


4.8% 



TABLE III: Dark energy baseline constraints and degradations for the fiducial model (tr^gd = 0.03(1 + z), N z = 20 and Sz p = 0.1) 
and deviations from the fiducial photo-z parameters. Baseline constraints are for <r(z blas ) = o-(cr^) = and relative degradations are for 
CT (z bias ) = 0.003,/ N z /20 and u{a 2 z ) = (0.03) 2 ^N z /20. The priors are scaled to reflect constraints per fixed bin Sz = 0.1 of true redshift. 



tice from Tabic IIIII that the baseline constraints do not 
change much by doing this changes in the fiducial model 
but the degradations relative to the baseline are affected 
by a series of competing effects. 

Decreasing the number of photo-z parameters to N z = 
10 while keeping the photo-z bin size at 8zf = 0.1 (i) 
partially decreases some degeneracies between dark en- 
ergy and the fewer photo-z parameters and (ii) makes the 
photo-z parameters smoother functions of redshift, mak- 
ing photo-z perturbations more nonlocal (less restricted 
to a particular photo-z bin) . The latter effect causes (a) 
the counts sensitivity to z b,as to increase because it de- 
pends on the value of z blas at the bin borders, and (b) the 
sensitivity to a z decreases because the number of objects 
scattered up, controlled by the values of a z at the bin bor- 
ders, more efficiently compensates the number of objects 
scattered down, controlled by a z at the bin center. In 
the case of counts only and perfect mass the dark energy 
degradations are weaker than the fiducial case whereas 
in the case of unknown masses and self-calibrations they 
are stronger. 

Increasing the bin size to 8z p = 0.2 but keeping 
N z = 20 (i) produces a coarser probe of cluster mass 
function but (ii) reduces the Poisson noise in the counts, 
therefore increasing the importance of sample variance 
and clustering self-calibration. Notice though that in all 
cases absolute errors increase with dz p since subdividing 
the data can only help constraints. 

If in addition to using coarser bins of 5z p = 0.2 we 
simultaneously decrease the number of photo-z parame- 
ters to N z = 10 so as to still have one photo-z param- 
eter per photo-z bin, the degradations in all cases be- 
come slightly better than in the fiducial model, although 
around slightly worse baselines. That shows that our re- 
sults arc relatively robust to the particular binning choice 
if we keep the number of photo-z parameters per bin of 
true redshift constant. In this case, coarser photo-z bins 
lose shape information on the redshift distribution but al- 



lows the fewer photo-z parameters to be better known by 
the same factor that they were reduced in number. Since 
the best achievable constraints are worse with such wider 
bins and the power of cluster counts surveys lies exactly 
in their ability to split the counts in a reasonable num- 
ber of redshift and mass bins, photo-z bins much coarser 
than these are unlikely to be adopted in future cluster 
survey analyses if they are to be competitive dark energy 
probes. 



E. Training Set Requirements 

The photo-z requirements for a fixed dark energy 
degradation ultimately translate into requirements on the 
size of the calibration training set (with known rcdshifts), 
from which we infer the photo-z error parameters. With 
the assumption made here that the cluster photo-z er- 
rors have a Gaussian distribution, the uncertainties in 
the photo-z parameters z hlas (zi) and cr^(zi) in a given 
redshift bin i are related to the number N spec (zi) of spec- 
troscopic training set calibrators in that bin by 

a(z bias (z 4 )) = - 1 (30) 

y N s p ec {Zi) 

a(a 2 z ( Zi )) = o*{zi)J 2 ■ (31) 
y N spec \Zi) 

We are interested in the values of N spcc necessary in 
order not to degrade dark energy constraints appreciably. 
There are basically two ways to go about this question, 
which depend on the choice of the spectroscopic calibra- 
tors. The first case corresponds to choosing the calibra- 
tors to be clusters with at least one member galaxy whose 
spectroscopic redshift has been measured, in which case 
N spcc corresponds to the number of such spectroscopic 
clusters and <r z (zi) is the cluster photo-z scatter in the 
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rcdshift bin. This choice is natural, but may suffer from 
the fact that a small number of such clusters are likely 
to be available in a given survey. The second choice then 
tries to compensate for that by including more objects 
in the calibration process. Since most clusters contain 
large numbers of red elliptical galaxies and optical cluster 
finder methods are typically tuned to find exactly such 
clusters, one can try to use not only the cluster spec- 
troscopic redshifts, but the redshifts of all spectroscopic 
elliptical galaxies including field galaxies. Although this 
would make use of a much larger sample of objects, they 
would typically have a larger photo-z scatter and the as- 
sumption that all red galaxies (in the field and in clusters) 
share the same photo-z error properties would have to be 
checked. In any case, the first choice is general enough 
in the sense that we can think of the red galaxy sample 
as being used, not directly in the cluster photo-z param- 
eters calibration, but in the estimation of these parame- 
ters. From that perspective, having a large spectroscopic 
sample of red galaxies helps calibrate the photo-z's of 
cluster member galaxies and decrease the cluster photo- 
z scatter, since each cluster member provides roughly an 
independent photo-z measurement of the cluster photo-z. 
Even though some cluster finding methods, such as those 
based on the SZ signal and weak lensing shear, do not 
preferentially detect clusters with large numbers of red 
galaxies, most future cluster surveys will have a spectro- 
scopic follow-up of a sample of the detected clusters. For 
that reason, in the following analysis we will take the first 
choice of rcdshift calibrators, i.e. spectroscopic clusters. 

Recall in Fig. [2] we considered fixed degradations in 
the dark energy equation of state w, given priors <7(,z blas ) 
and ct(ct^) which were constant in redshift. In that case, 
each value of d w requires er(z blas ) and <j{a^) to be lower 
than certain values, which in turn produce lower limits on 
N sp ec( z i)- For illustrative purposes let us fix d w = 10% 
and consider the fiducial model with counts information 
only, in which case Fig. [2] requires cr(z blas ) < 0.003 and 
y/a(o%) < 0.03. The condition on cr(2; blas ), which is more 
stringent than the one on a{a 2 z \ requires the number of 
calibrators in bins of Szi = 0.1 to be N spcc (zi = 0.5) > 
110 and N spec (zi = 1.95) > 870. However, all these 
numbers are excessively high due to our requirement of 
constant a(z hlas ) and a(a z ) in each bin, including the 
high redshift ones that have little impact on the dark 
energy constraints. 

To remedy this problem, let us instead take N spcc to 
be a fixed fraction of the total number in a given bin 



-^spcc(^) 



(32) 



We can then study the requirements on dark energy 
degradations as a function of a. We should keep in mind 
though that, obscrvationally, it is even harder to obtain 
spectra for the few clusters at high redshifts compared 
to those (also few ones) at very low redshift, because 
at higher redshifts the 4000A break of elliptical galaxies 
leaves the optical filter coverage. Therefore in practice a 
will probably be a decreasing function of rcdshift. 
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FIG. 8: Dark energy parameters degradation in the dark energy 
equation of state as a function of the fraction a of clusters with 
spectroscopic redshifts for the fiducial model, WMAP3 cosmology 
and photo-z scatter of Tzlfid = 0.02, using counts information only 
and self-calibration. 



In Fig.[8j we show the degradation d w as a function of a 
in the fiducial model and also for the WMAP3 cosmology 
as well as the case with (t z \sa = 0.02. Notice that even 
in the limit where all clusters have photo-z's (a = 1), be- 
cause the number of these clusters is finite, we still have 
some degradation compared to the case where the photo- 
z parameters have zero uncertainties. Even though this 
intrinsic degradation can be as high as ~ 10%, as shown 
in the top panel of Fig. [8j if a was actually equal to 1, 
one would not use photo-z's in the first place, since all ob- 
jects would have perfect rcdshift measurements. There- 
fore actual degradations are smaller than those shown in 
Fig. [5] for very high values of a. Also, in practice only 
values of a that realize integer numbers of clusters are 
allowed. However, since we are considering a constant 
in redshift, we used continuous values of a, even those 
that give iV S pec < 1. Thus, actual degradations for very 
low values of a are larger than those shown in Fig. [8j 
especially for WMAP3. 

Ignoring these caveats just mentioned, we see from 
the top panel in Fig. [8J that the WMAP3 case has a 
smaller degradation than the fiducial model when em- 
ploying counts only, consistently with the results of 
ilVD 11 Moreover reducing the photo-z scatter has an 
important impact on the w degradation as expected. 

From the bottom panel of Fig. [5J we see that self- 
calibration tends to make the training size requirements 
smaller in most cases. Interestingly, d w in both the fidu- 
cial model and the WMAP3 cosmology are nearly iden- 
tical as a function of a. That happens because for fixed 
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values of a(z hlas (zi)) and a(cr^(zi)), whereas w degra- 
dations are typically weaker for WMAP3, the numbers 
m, of existing clusters are smaller and consequently the 
number of spectroscopic calibrators are reduced by the 
same factor, making it harder to achieve these photo-z 
accuracies. 

These results highlight the importance of having the 
best photo-z estimators available, which allow for the 
photo-z scatter to be as small as possible; this has a di- 
rect impact in reducing the required calibration set size. 
In the context of self-calibration, having spectroscopic 
redshifts for only ~ 1% of the clusters is enough to keep 
d w Hz 30%, whereas in the case of counts only, this spec- 
troscopic sample needs to contain at least ~ 8% of the 
clusters. For a more stringent requirement of d w < 20%, 
with self-calibration it is necessary to have spectroscopic 
redshifts for ~ 4% of the clusters. Quantitatively, this 
requires ~ 1300 spectroscopic clusters in the fiducial 
model and ~ 300 in the WMAP3 cosmology. Given that 
the training set size requirements for other cosmological 
probes such as weak lensing and baryon acoustic oscil- 
lations arc on the order of hundreds of thousands, only 
a few percent of these galaxies need to be in clusters in 
order for those training sets to serve also as cluster cali- 
bration samples, without additional observational effort. 
Complementing the existing training sets with spectro- 
scopic follow-up observations would then further improve 
dark energy constraints and allow for degradations at the 
percent level. 



F. Caveats 

There are a number of caveats associated with our re- 
sults. We have taken a constant M°^ s whereas in reality 
we expect it to be cosmology and redshift dependent for 
a SZ survey like SPT (although less than in an X-ray 
survey [H|). Our constant fiducial M$* can be thought 
of as an effective mass threshold that produces the same 
number counts as the true threshold. If the effective mass 
threshold in a real survey deviates considerably from the 
one considered here, the baseline dark energy constraints 
will change, and so will the degradations with respect to 
this baseline. As discussed in fjV] a change in MP^ has 
an effect similar to a change in ag or in the survey area, 
namely the change in the total cluster counts. 

The variation of the observable-mass relation with cos- 
mology may also open up parameter degeneracies. For 
our fiducial model we took the observable-mass distri- 
bution to be parametrized by lnM blas and crf nM with 
a fixed functional form in redshift. If this form instead 
mimicked the effect of a cosmological parameter, a de- 
generacy would appear. A complete treatment of these 
issues, which is beyond the scope of this work, would 
require the choice of a specific observable and extensive 
simulations to obtain the true parameter values and their 
distributions. 

The potential impact of degeneracies can be assessed 



by considering the case where In M blas and of n M are al- 
lowed to be arbitrary functions of redshift ( W Dp . In this 
case the baseline dark energy constraints degrade by up 
to 80% relative to the fiducial model. Degradations on 
top of that due to uncertainties in z blas and o\ become 
weaker and stronger respectively. 

If one further includes an arbitrary mass dependence 
of these parameters, the techniques explored here are 
not sufficient to provide a good degree of self-calibration, 
even with perfect redshifts. Our results would still be 
valid as long as trends in mass are well-known close to 
the mass threshold Q • Recent simulations [13, HH, H3 
suggest that observable-mass relations can actually be 
parametrized with scaling relations involving few param- 
eters as assumed here. Since our main purpose here is 
to illustrate the relative degradation of dark energy con- 
straints due to photo-z errors, the range of models con- 
sidered here should suffice. 

Finally, we assumed the photo-z error parameters to 
follow a Gaussian distribution and depend on redshift 
but not mass. This assumption does not hold for field 
galaxies binned only in redshift due to the mixing of var- 
ious galaxy types. Even though we expect Gaussian er- 
rors to be a good approximation for clusters, which have 
mostly red elliptical galaxies and lower photo-z scatter, 
possible non-Gaussian features induced by membership 
contamination may require proper modeling. 

The photo-z parameters, in particular the photo-z scat- 
ter is expected to have some mass dependence. More 
massive clusters have more galaxies to average over when 
computing the cluster photo-z, allowing the photo-z scat- 
ter to decrease. We envision that the redshift/mass de- 
pendence of photo-z parameters can be modeled from 
mock catalogs and cast in simple scaling relations, likely 
depending on the particular cluster finder/survey. In the 
lack of any specific modeling, we chose to keep the red- 
shift dependence general and ignore the mass dependence 
of the photo-z parameters. Lastly, we did not take into 
account possible errors in the theoretical cluster-mass 
function, which might introduce parameter biases or in- 
teract differently with photo-z errors. 



VI. DISCUSSION 

The use of cluster abundance as a competitive cos- 
mological probe requires strict knowledge of both the 
cluster observable-mass relations and photo-z's. Self- 
calibration techniques help break degeneracies between 
cosmology and cluster masses and, at the same time, de- 
crease redshift knowledge requirements. When external 
mass calibrations are available, self-calibration may be- 
come even more important, assisting monitor redshifts 
and providing interesting consistency checks. Conversely, 
good photo-z knowledge is required in order to extract 
optimal constraints from self-calibration methods. 

We studied the effect of photo-z uncertainties in dark 
energy constraints assuming a fiducial model with sur- 
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vey specifications similar to the SPT, a constant mass 
threshold M° h bs , a mass bias lnM blas parametrized by a 
power law in redshift, a mass variance of nM with a cubic 
redshift evolution and photo-z parameters assumed to be 
general functions of redshift but not mass. 

For this fiducial model it will be necessary to have 
priors of cr(z bias ) < 0.003 and y r a\af} < 0.03 in the 
case of perfect masses and counts information in order 
not to degrade dark energy constraints by more than 
*~ 10%. In the case of self-calibration of the observable- 
mass relation, these priors uncer tainties should be kept 
at er(z bias ) < 0.004 and y/aipfj < 0.04. These require- 
ments become weaker/stronger if we increase/decrease 
the number of nuisance parameters describing masses and 
redshifts or if we consider less/more ambitious surveys or 
cosmologies yielding less/more clusters. 

In order to achieve these requirements, it is impor- 
tant not only to use the best general photo-z techniques, 
but also to understand the cluster finding selection func- 
tions and contamination fractions; otherwise even if the 
photo-z methods work well for field galaxies, the clus- 
ter redshifts and masses might be misestimated. Large 
and representative training sets help decrease the red- 
shift uncertainties and keep the dark energy constraints 
close to their baseline values and also help self-calibration 
be more effective. From the self-calibration perspective, 
simulations play a very important role. Self-calibration 
techniques are more effective when we can parametrize 
the observable-mass relation with a small number of pa- 
rameters, typically with simple scaling relations. Simula- 
tions do not need to determine exact values of these pa- 
rameters, but only provide a good confidence that these 
simple relations are stable to theoretical uncertainties 
and have relatively small theoretical scatter. If that is 
the case one can extract these parameter values from 
self-calibration methods, which use information that nat- 
urally comes in surveys without requiring additional ob- 
servational efforts. 

Photo-z methods are expected to be ever improv- 
ing during the evolution of ongoing and future surveys 
such as the Panoramic Survey Telescope and Rapid Re- 
sponse System (PanSTARRS) [32| , the SPT [29[, the 
Dark Energy Survey (DES) [331, the Supernova Acceler- 
ation Probe (SNAPJ_J34[ and the Large Synoptic Survey 
Telescope (LSST) J35(. There are good prospects for the 
existence of large and representative training sets that 
will allow these surveys to be reasonably well calibrated. 
Likewise we expect simulations to improve our knowl- 
edge of the dark matter halo mass-functions and of the 
various observable-mass relations. External calibrations 
from weak lensing are likely to also play an important role 
in mass determinations. Finally, self-calibration methods 
can further help reduce requirements and provide impor- 
tant consistency checks that theoretical assumptions have 
been or not realized in the observations. All these aspects 
will be important and complementary in consistent cos- 
mological analyses of future cluster surveys. 
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APPENDIX A: WINDOW FUNCTION 

Recall we used spherical coordinates x = (r, 9, <f>) 
to parametrize the position vector describing the win- 
dow volume element such that d 3 x — r 2 dfldr with 
dfl = sm9ddd(j). The cell defines a solid angle given 
by Afl = tt9 2 . in the small angle approximation. For 
convenience, we will also use cylindrical coordinates to 
parametrize the same vector x = (r,p,(f>), where r(z) 
is still the angular diameter distance to redshift z and 
now the perpendicular coordinates (/?, (j)) parametrize the 
angular extension of the window. In these coordinates 
we have d 3 x ~ pdpd(f>dr and the cell defines a exten- 
sion of p s = r9 s . We split the total window function 
into Wi(x) = W] h i (p,(j))F i (z). Here and below the su- 
perscript "th" denotes a top hat function whose value is 
unity within the pixel and zero outside of the pixel. 

The volume-normalized Fourier transform of the win- 
dow is given by 



3 xe lkx 



VF,(x) 



-ini(k±)W^(fc||), 



where the window volume is 



Vi 



and we defined 



d 3 xWi{x.) = An / drr 2 Fi(z) 



(Al) 



(A2) 



W£(kj.) ee j ' dpdtpe^^+WfiM), (A3) 
W ¥ (k() ee f dre ik « r F t (z). (A4) 



Note that the normalized Fourier window Wj(k) — ► 1 as 
k^0. 

Similarly, in the case of a top hat total window we 
would have 



1 



w^(k) = _wS(k J .)wff(ftn) 



with 



V? 



th _ / j3^,TT/th/ 



(A5) 



d 3 xW^(x) = Ail I drr 2 W^(r) . (A6) 



From Eqs. (|A1|) and (|A5[) we may write 



(A7) 
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and we can rewrite Eq. (|A4[) as 

dz 



In the limit where r(z) and the photo-z parameters 
and a\ do not change appreciably inside the photo-z 



v bias 



/ d z c ik u r(z) p f z \ (A8) ^^ n ' wm( ^ ow volumes in the absence (V/ h ) and pres 

H(z) ence of photo-z errors (VA roughly coincide V, « V, th o 



We will assume that H(z) does not change apprecia- 
bly in the photo-z bin and, because Fi(z) quickly drops 
outside the bin, the values of r that contribute to the 
integral are restricted to the range r = r{ ± 5ri/2 with 
~ $ z i ' /Hi and Hi = H(zf), so that 

wdh) * Wi j dzeik[lZ/HiF ^) = ■ ( A9 ) 

where we defined k z ,i = kn/H. Note that Fi(z) can be 
expressed as a convolution of the Gaussian selection with 
a top hat window in rcdshift 

Fi(z)= [ dz phot Wt(z phot )p(z phot \z) . (A10) 



The convolution theorem gives Fi(k z ,i) 
W^(k zi )p(k zi ) so that Eq. §M§ becomes 



W^(k^)p(k Zt i) 



H 



(All) 



By a similar argument that led to Eq. (|A9[) we have that 
Wl h (fc 2i4 ) ~ HiW$(k\\) and Eq. ((13) becomes 



Wi(k) V? h 

= —p(k z ,i) 



W, th (k) v 



(A12) 



ence of photo-z errors (Vi) roughly coincide V, « V/ 
rfSri, i.e. the photo-z errors distort the volume element 
but do not change its value. Therefore we have 

Wi(k) = W} h (k)p(k Zti ). (A13) 

The Fourier transform of the top hat window is [2(| 

ifcnn sin(fc||<5r J /2) Ji(fc_Lr« 



M/* h (k) = 2e J 

fc||5ri/2 fc±ri6l s 

and for the Gaussian distribution we have 

p(k zA ) - e *.,.»Y" e -<rt,J*. 

Combining these results, we obtain 

„bias 



(A14) 



(A15) 



Wi(k) = 2exp 



ikn 



H 



exp 



~2Hf 



sin(fc| | (5r i /2) J^jk^rA) 



for the window function. 
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